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Abstract: Using an 18-year dataset of arrival dates of 65 species of Maine migratory 
breeding birds, I take a deeper view of the data to ask questions about the shapes of the 
distribution. For each year, most species show a consistent right-skewed pattern of distribution, 
suggesting that selection is stronger against individuals that arrive too early compared to those 
that arrive later. Distributions are consistently leptokurtic, indicating a narrow window of 
optimal arrival dates. Species that arrive earlier in the spring show higher skewness and 
kurtosis values. Nectarivorous species showed more pronounced skewness. Wintering area 
did not explain patterns of skewness or kurtosis. Deviations from average temperatures and 
the North Atlantic Oscillation index explained little variation in skewness and kurtosis. 
When arrival date distributions are broken down into different medians (e.g., 5% median 
and 75% median), stronger correlations emerge for portions of the distribution that are 
adjacent, suggesting species fine-tune the progress of their migration. Interspecific 
correlations for birds arriving around the same time are stronger for earliest migrants 
(the 25% median) compared to the true median and the 75% median. 

Keywords: birds; migration; skewness; kurtosis; arrival dates; weather; climate change; 
phenology 



1. Introduction 

Abundant evidence for global climate change exists in monotonic increases in carbon dioxide 
concentrations, melting of polar ice caps with concomitant sea level rise, and record hot temperatures 
around the world [1^]. Pheno logical data add to this body of evidence [5-13]. Changing migration 
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schedules of birds provide some of the most compelling indications that these climatic changes are 
influencing the natural world. Birds are unique barometers of climate change because of the rich 
database on migration contributed in large part by lay citizens [5,6,12,13]. Recent reviews of the large 
quantity of data leave little doubt of a fundamental impact of climate change on bird migration [14-16] 
although some of the specific mechanisms of migratory control and behavior have been recently 
challenged [17]. 

Over the past 18 years, I have been coordinating a citizen-science effort to document the first arrival 
dates of Maine migratory breeding birds across the state [18-20]. Previous comparisons of these 
modern data to data taken over a century ago from the Journal of the Maine Ornithological Society 
(published from 1899 until 1911) provide mixed evidence of earUer contemporary arrivals for some 
but not all species [21]. Over the eighteen-year period of the study, arrival dates are associated with 
fluctuations in mean springtime temperature as well as by the strength of the North Atlantic Oscillation 
(NAO) [19,20]. 

In this contribution, I take a longer view of the data. Rather than concentrating on the impact of 
various climatic variables on arrival dates, I explore the statistical distribution of arrival dates reported 
by all observers. As argued by [19], migratory birds are under competing pressures that result in 
stabilizing selection. On the one hand, birds, especially males, are being pushed to arrive as early as 
possible on the breeding grounds to stake out the most favorable breeding territories. On the other 
hand, birds are pressed to delay arrival in Maine until sufficient food is available to support them. We 
therefore expect that the distribution of arrival dates might yield a normal distribution. However, a 
positive or negative skew to the data could indicate that one of the two selective pressures is stronger 
than the other. I examine the distributions of arrival dates for each species and year for skewness and 
kurtosis. Kurtosis, a measure of the peakedness of a distribution, might be influenced by these 
selective pressures as well. A positive kurtosis (leptokurtosis) would indicate a prominent peak with 
particularly narrow tails, indicative of strong stabilizing selection. I ask if skewness and kurtosis are 
related to variation in springtime temperature, to the NAO, to foraging types of the birds, and to 
wintering area of the birds. 

I also break down the distributions to see if the dates of arrivals of the earliest arriving individuals 
of a species can predict the arrival of the later arriving birds. I examine the interannual correlations of 
arrivals for 65 species of birds. With this deeper view, I am less focused on the responses of particular 
birds and more focused on discovering patterns by examining elements of the distributions of 
arrival dates. 

2. Materials and Methods 

The data for this study come from the 18 years (1994—2011) of data on arrival dates of Maine 
migratory breeding birds [18-20]. Each volunteer records the first arrival of over 100 species of 
migratory breeding bird and the Biophysical Region (Figure 1) in which each observation was 
made; [22] delineated these Biophysical Regions based on climate and vegetation patterns. Over 
44,000 records have been reported by over 200 volunteers across the state. Although I have data on 
145 species of migratory breeding birds, I restrict this study to an analysis of each of the 65 species 
that have a minimum of 20 arrival dates reported for each year of the study. For each spring, I 
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calculated the departure of mean temperature (TempDep) from the 100-year mean. I obtained the 
100-year mean by averaging data from 15 different meteorological stations across the state. 
I determined TempDep for March, April and May. The appropriate TempDep value was used for 
species that mostly arrive within a particular month. However, some species have arrival dates that 
span two months (e.g., late March into early April for Killdeer). For such species that arrive in late 
March into early April, I averaged the March and April TempDep values. Similarly, I averaged April 
and May for those species that arrive in late April into early May. I also used the NAO index for each 
year of the study. 

Figure 1. Map of Biophysical Regions of Maine, based on [22]. Map from the Maine 
Office of Geographic Information Systems (MEGIS). 




Birds were categorized into one of nine foraging types (Table 1). Most of the foraging types are 
straightforward [23]. Gleaners take insects from leaves, granivores take seeds (mostly from the 
ground) while ground predators take animal prey from the ground. Scansorial predators take insects 
and other arthropods from the bark of trees. 
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Table 1. Foraging types, skewness, kurtosis and major migration period of the 65 species 
used in this study. The skewness and kurtosis values are the mean of the 18 yearly values. 
An asterisk after particular skewness or kurtosis values indicates that the value is 
significantly different from zero in Wilcoxon rank tests using the Bonferroni correction to 
insure an experiment-wise /7-value of 0.05. 



Species 


Foragiiig Type 


Skewness 


ICurtosis 


Migration 
Period 


Ring-necked Duck (Aythya collaris (Donovan)) 


Aquatic 


0.240 


3.354 * 


April 


Great Blue Heron {Ardea herodias (Rackett)) 


Aquatic 


0.620 * 


4.667 * 


April 


Turkey Vulture (Cathartes aura (L.)) 


Scavenger 


0.509 * 


3.757 * 


March/ April 


Osprey {Pandion haliaetus (L.)) 


Raptor 


0.699 * 


4.687 * 


April 


Broad-winged Hawk (Buteo platypterus Viellot)) 


Raptor 


-0.390 


4.249 * 


April 


American Kestrel {Falco sparverius L.) 


Raptor 


0.449 * 


4.001 * 


April 


Killdeer {Charadrius vociferus L.) 


Ground predator 


0.873 * 


4.392 * 


March/ April 


Spotted Sandpiper (Aclilis iiiaciilariiis (L.)) 


Ground predator 


0.161 


2.965 * 


May 


Wilson's Snipe {Gallinago delicata Ord) 


Grotmd predator 


0.476 * 


2.861 * 


April 


American Woodcock (Scoloplax minor Gmelin) 


Grotmd predator 


1.231 * 


5.490 * 


March/ April 


Chimney Swift {Chaetura pelagica (L.)) 


Aerial 


0.523 


3.863 * 


May 


Ruby-throated Hummingbird {Archilochus colubris (L.)) 


Nectarivore 


0.108 


6.901 * 


May 


Belted Kingfisher (Megaceryle torquata (L.)) 


Aquatic 


-0.636 


5.357 * 


April 


Yellow-bellied Sapsucker (Sphyrapiciis variiis (L.)) 


Scansorial 


1.479* 


3.578 * 


April 


Northern Flieker {Colaples aiiraliis (L.)) 


Seansorial 


-0.783 


7.253 * 


April 


Eastern Wood-Pewee (C. virens (L.)) 


Aerial 


-0.729 


6.373 * 


May 


Least Flycatcher (E. minimus (Baird and Baird)) 


Aerial 


0.283 


3.745 * 


May 


Eastern Phoebe (Sayornis phoebe (Latham)) 


Aerial 


1.907 * 


5.463 * 


March/ April 


Great Crested Flycatcher {Myiarchus crinitus (L.)) 


Aerial 


0.582 


4.525 * 


May 


Eastern Kingbird {Tyr annus tyrannus (L.)) 


Aerial 


-0.518 


6.045 * 


May 


Blue-headed Vireo {Vireo solitarius (Wilson)) 


Leaf-gleaner 


0.478 


3.832 * 


April/May 


Warbling Vireo {V. gilvus Viellot) 


Leaf-gleaner 


0.533 * 


2.883 * 


May 


Red-eyed Vireo (V. olivaceus (L.)) 


Leaf-gleaner 


0.335 


5.247 * 


May 


Tree Swallow (Tachycineta bicolor (Viellot) 


Aerial 


0.945 * 


4.166* 


April 


Bam Swallow {Hirundo rustica L.) 


Aerial 


0.173 


3.037 * 


May 


House Wren {Troglodytes aedon Viellot) 


Ground predator 


0.277 


3.239 * 


May 


Winter Wren (Troglodytes troglodytes L.) 


Ground predator 


0.627* 


2.792* 


April 


Ruby-crowned Kinglet (Regulus calendula (L.)) 


Leaf-gleaner 


0.454 


3.725 * 


April/May 


Eastern Bluebird (Sialia sialis (L.)) 


Grotmd predator 


0.225 


3.708 * 


April 


Veery (Catharus fuscescens (Stephens)) 


Grotmd predator 


-0.181 


3.808 * 


May 


Hermit Thrush (C. guttatus (Pallas)) 


Grotmd predator 


0.301 


4.693 * 


May 


Wood Thrush {Hylocichla miistelina (Gmelin)) 


Ground predator 


0.175 


3.891 * 


May 


Gray Catbird (Dumetella carolinensis (L.)) 


Ground predator 


-0.637 


8.439 * 


May 


Brown Thrasher (Toxostoma nifiim (L.)) 


Ground predator 


0.217 


3.579 * 


May 


Ovenbird (Seiurus aurocapilla (L.)) 


Leaf-gleaner 


0.216 


5.452 * 


May 


Northern Waterthrush (Parkesia novaeboracensis (Gmelin)) 


Leaf-gleaner 


0.367 


3.585 * 


May 


Black-and-white Warbler {Mniotilta varia (L.)) 


Leaf-gleaner 


0.653 


3.838 * 


April/May 


Nashville Warbler (Oreothlypis ruficapilla (Wilson)) 


Leaf-gleaner 


0.804* 


4.546 * 


May 
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Table 1. Cont. 



Species 


Foraging Type 


Skewness 


Kurtosis 


Migration 
Period 


Common Yellowthroat (G. tncnas (L.)) 


Leaf-gleaner 


0.928 * 


5.056 * 


May 


American Redstart (Setophaga ruticilla (L.)) 


Leaf-gleaner 


0.448 


5.434 * 


May 


Northern Parula (5. amencana (L.)) 


Leaf-gleaner 


0.098 


5.056 * 


May 


Magnolia Warbler (5. magnolia (Wilson)) 


Leaf-gleaner 


0.262 


3.562 * 


May 


Blackburman Warbler (S. jusca (Forster)) 


Leaf-gleaner 


0.628 


4.514 * 


May 


Yellow Warbler {S. petechia (L.)) 


Leaf-gleaner 


0.470 


4.187 * 


May 


Chestnut-sided Warbler (S. pensylvanica (L.)) 


Leaf-gleaner 


0.629 


4.823 * 


May 


■ -V 1 1^1 ^1 1 "\ 1 11711 /f~T / / 1 ■ W 

Black-throated Blue Warbler (S. caerulescens (Gmelm)) 


Leaf-gleaner 


0.652 * 


3.838 * 


May 


Palm Warbler {S. palmarum (Wilson)) 


Leaf-gleaner 


0.579 


3.824 * 


April 


• TTT 11 • /TT T"*1 

Pine Warbler(5. pinus (Wilson)) 


Leaf-gleaner 


0.651 


4.208 * 


April 


T 7" 1 1 1 TT T 1 1 / A /T "W 

Yellow-rumped Warbler (S. coronata (L.)) 


Leaf-gleaner 


-0.835 


8.506 * 


April/May 


■ -Y 1 1.1 .1 117 11 /" o ■ 1 ■ \\ 

Black-throated Green Warbler {S. virens (Gmelm)) 


Leaf-gleaner 


0.667 * 


4.061 * 


May 


l 117 11 / 111- 1 - /TW 

Canada Warbler {Caraellina canadensis (L.)) 


Leaf-gleaner 


0.320 


3.429 * 


May 


Chipping Sparrow (Spizella passerina (Bechstein)) 


Granivore 


-0.254 


4.446 * 


April 


Savannah Sparrow (Fasserculus sandwicnensis (Gmelin)) 


Granivore 


0.014 


3.676 * 


April 


Song Sparrow (Melospiza melodia (Wilson)) 


Granivore 


-0.792 


8.560 * 


Mar/April 


Swamp Sparrow (M georgiana (Latham)) 


Granivore 


0.094 


3.426 * 


April 


TT T1 'j j1 j 1 n / J • 7 • TT • TT* 1* 

White-throated Sparrow {Zonotrichia albicollis (Gmelin)) 


Granivore 


-1.986 * 


7.842 * 


April 


Scarlet Tanager {Piranga olivacea (Gmelin)) 


Leaf-gleaner 


0.651 * 


3.746 * 


May 


Rose-breasted Grosbeak (Pheucticus ludovicicmus (L.)) 


Granivore 


0. /o4 


c Am * 
b.\r)l ^ 


May 


Indigo Bunting {Passerina cyanea (L.)) 


Granivore 


0.157 


3.026 * 


May 


Bobolink {Dolichonyx oryzivorus (L.)) 


Ground predator 


0.605 


7.210* 


May 


Red-winged Blackbird (Aegelaius phoeniceus (Mtiller)) 


Ground predator 


1.483 * 


6.222 * 


Mar 


Eastern Meadowlark {Stumella magna (L.)) 


Ground predator 


1.992 


2.583 * 


April 


Common Grackle {Quisculus quiscala (L.)) 


Ground predator 


1.426 * 


7.179* 


Mar 


Baltimore Oriole {Icterus galbula (L.)) 


Leaf-gleaner 


0.010 


7.857 * 


May 



I classified each of the 65 species according to its major wintering area: South America, Central 
America, Caribbean, or southeastern United States. Central and South America along with Caribbean 
were combined into a Long category for the long-distance migration those birds undertake. Long contrasts 
with Short, referring to those birds wintering in the United States with a less arduous migration. 

For each species-year combination, I used Stata 12 statistical software to calculate the skewness and 
kurtosis [24]. Because skewness and kurtosis data are not normally distributed, non-parametric tests 
were used to test for differences among factors. I used the Wilcoxon signed-rank test to determine if 
the skewness and kurtosis for each species differed from a normal distribution by testing the observed 
values against a value of 0 (the value of skewness and kurtosis in a normal distribution). For both the 
skewness and kurtosis analyses, I used the Bonferroni correction to minimize Type I errors from 
performing multiple analyses [24]. The correction was set to ensure that an experiment -wise value of 
p = 0.05 was obtained by setting the critical /?-value for each test as 0.0007, obtained by dividing 0.05 
by 65 (the number of tests performed). Using the non-parametric Kruskall-Wallis rank test, I examined 
the relationship between arrival month, foraging type, migration distance and wintering area on both 
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kurtosis and skewness. I used the Bonferroni correction for these four tests, setting the critical value 
of each test io p = 0.0125, yielding an experiment -wise critical /?- value of 0.05. 1 used linear regression 
analysis to examine the ability of TempDep and NAO index to predict kurtosis and skewness. The 
critical />-value of 0.0007 was used here as well to produce an experiment -wise critical value of/? = 0.05. 

For the distribution of each species-year combination, I calculated five metrics: the 5% median, the 
10% median, the 25% median, the true median (50%) median) and the 75% median. I used Pearson 
product-moment correlations to evaluate the pair-wise relationship between each type of median across 
the 18 years of the study. 

Using the values of the 25% median, 50% median and 75% median, I used Pearson correlation tests 
to test for significant relationships between pairs of species that were migrating in the same time 
interval. Thus, I developed one table of all possible combinations of species arriving in a particular 
month. This analysis tests for the independence of migration schedules among the species in the study. 

3. Results 

The quantity of arrival data for each Biophysical Region (Figure 1) in this study refiects the human 
population density. The most observations came from the Central Interior (31.6% of all observations), 
the South Coastal Region (21.8%), the Midcoast Region (17.3%), the Penobscot Bay Region (8.2%) 
and the East Coastal Region (7.9%). Most of the remaining data (13.2%) were spread equally across 
the Southwest Interior, the Western Foothills, the Eastern Lowlands and the Eastern Interior. 

Table 1 presents the mean values of skewness and kurtosis for the 65 species examined in this 
paper. In a symmetrical distribution, skewness would be equal to zero. Fifty-four of the 65 species 
have a positive value of skewness. Using Wilcoxon signed-rank tests, 20 of the species show skewness 
values that are significantly different from zero. White-throated Sparrow is the only species that has a 
statistically significant negative skewness. Figure 2 presents the distribution of arrival dates of 
American Woodcock for 1994 as a typical example. This distribution has a positive skewness with the 
longer tail to the right. 

Figure 2. Distribution of first arrival dates for American Woodcock for 1994. The 
distribution is right-skewed with some leptokurtosis (more pronounced peak and thinner 
tails than seen in a normal distribution). 

American Woodcock 




80 90 100 110 

Day of Year 
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Wilcoxon signed-rank tests indicate every species has a significant positive kurtosis (a value of zero 
is expected for a normal distribution). The positive kurtosis is termed leptokurtosis, indicating a 
pronounced peak in the middle of the distribution with thinner than expected tails to the distribution. 

I used a Rruskal-Wallis rank test to determine if foraging type (Table 2) was associated with the degree 
of skewness. The differences were not significant (X =21.81,/' = 0.01). Scavengers and leaf-gleaners 
had the highest positive skewness and raptors had the lowest value. Similarly, a Kruskal-Wallis rank 
test, examining the effect of foraging type on kurtosis was also significant (X = 65.64, p < 0.01). 
Nectarivores were particularly leptokurtic with aquatic feeders and scavengers showing less 
pronounced leptokurtosis. 

Table 2. Relationship of foraging type on kurtosis. A Kruskal-Wallis rank test indicated 
significant differences among feeding types. 



Foraging type 


Skewness 


Kurtosis 


Aerial insectivore 


0.38 


4.65 


Aquatic feeder 


0.35 


3.09 


Gleaner 


0.45 


4.59 


Granivore 


0.32 


5.48 


Ground predator 


0.31 


4.41 


Nectarivore 


0.06 


7.38 


Piscivore 


0.22 


4.13 


Raptor 


0.03 


4.90 


Scansorial predator 


0.35 


5.42 


Scavenger 


0.51 


3.76 



A Kruskal-Wallis rank test was used to test the statistical effect of migration month on skewness 
and kurtosis (Table 3). A significant effect was found for skewness (X^ = 49.70, p < 0.01) and kurtosis 
(X^ = 30.43, p < 0.01). Early arriving species (Red-winged Blackbird and Common Crackle) 
were associated with the highest skewness and highest kurtosis. The May arrivals showed a slight 
negative skew. 

Table 3. The relationship of migration period on skewness and kurtosis. Means are 
provided. Migration period had significant effects on both measurements in Kruskal-Wallis 
rank tests. 



Migratioii period 


Skewness 


Kurtosis 


March 


1.45 


6.73 


March/April 


0.57 


4.65 


April 


0.26 


4.38 


April/May 


0.24 


4.57 


May 


-0.29 


4.83 



Kruskal-Wallis rank tests failed to show an effect of wintering area (and hence migratory distance) 
on skewness (X = 7.34, p = 0.02) and on kurtosis (X = 3.32, p = 0.65). Combining the three tropical 
regions into one region to contrast with the US region fails to show a significant effect as well 

2 2 

(X = 0.37,/? = 0.55 for skewness and X = 0.21, /» = 0.35 for kurtosis). 
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To examine the statistical effect of TempDep on skewness, I used regression analyses. With the 
Bonferroni correction, no significant relationship was found any of the 65 species (all values > 0.0007). 
Similarly, regression analyses were used to test the statistical effect of TempDep on kurtosis. None of 
the 65 regressions was significant (p > 0.0007 in all cases). 

No significant relationship of NAO index with skewness was detected by regression analysis for all 
65 species. Regression analyses also failed to find a relationship between NAO index and kurtosis. 
None of these 130 regressions jaelded a p-value below the Bonferroni-corrected value ofp = 0.0007. 

Table 4 provides a summary of the correlation of different portions of the distribution (5% median, 
10% median, 25% median, 50% median and 75% median). Percentages represent the fi-action of 
species showing a significant relationship. For instance, for the 10 March/ April migrants, the 5% 
median and 10% median are significantly correlated for all 10 species but only five of the 5% median 
to 25% median correlations are significant. A consistent pattern throughout is that adjacent medians 
are much more likely to be correlated than medians fi-om different parts of the distribution (5% median 
correlated with the 75% median, for example). 

Table 4. Percentage of species whose various medians show a significant Pearson 
correlation (p < 0.05). 



March- April (10 species) 





5% 


10% 


25% 


50% 


10% 


100 








25% 


50 


80 






50% 


20 


50 


10 




75% 


20 


30 


40 


100 



April (16 species) 





5% 


10% 


25% 


50% 


10% 


88 








25% 


38 


81 






50% 


25% 


43 


88 




75% 


25 


31 


44 


69 



April-May (6 species) 





5% 


10% 


25% 


50% 


10% 


83 








25% 


50 


67 






50% 


17 


17 


83 




75% 


0 


0 


17 


67 



May (3 1 species) 





5% 


10% 


25% 


50% 


10% 


87 








25% 


58 


84 






50% 


26 


55 


77 




75% 


6 


16 


32 


94 
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Figure 3 presents a summary of the interspecific correlation analyses. The 25% median, 50% 
median and 75% median were analyzed separately. Each species migrating in a particular month was 
analyzed with every other species for that month. A consistent pattern is that the percentage of 
significant correlations for the 75% median is always the lowest. For three of the four time intervals, 
the 25% median had the highest percentage of significant correlations. For the April-May interval, 
only six species were analyzed making general patterns harder to discern. 

Figure 3. Box plots showing the percentage of interspecific significant correlations 
ip < 0.05) using the 25% median, 50% median and 75% median. Data for March are not 
presented because only two species arrive then (Red-winged Blackbird and Common 
Crackle). The box plot for April-May is based on only six species, accounting for its 
collapsed appearance. 
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4. Discussion 

To seek effects of climate change on bird migration, ornithologists have used two sorts of data. A 
few studies have been conducted where the arrival of all members of a local population of a few 
species are recorded through a banding program or intensive censusing of territorial males [25,26]. 
This total population approach provides ideal data but the effort required is enormous. Instead, most of 
the data collected to examine the impact of global climate change on birds come from first arrival 
dates [13-15,17-21]. Such data are easily obtained and the contribution of amateurs to this database is 
significant. The disadvantage is that we are only sampling the left-hand portion of the distribution, the 
vanguard of the migratory wave [17]. 

The second approach was used in this study. The analysis in this contribution was designed to seek 
patterns of deviation from the normal distribution in the shape of the arrival date curves. Studies using 



Biology 2013, 2 



751 



arrival dates of migratory birds to seek evidence of global climate change have concentrated on the 
mean or median dates [2,5,6,8,9,1 1]. I argue that the shape of the distribution of arrival dates may offer 
additional insight. Patterns of skewness and kurtosis have been profitably quantified in studies using 
the total population approach [25,26]. 

The first result from this analysis is that the distributions of arrival dates have third (skewness) and 
fourth (kurtosis) moments that differ from a normal distribution. In the majority of cases, the 
distribution shows a right skew (Figure 2, Table 1). The distributions also show leptokurtosis, a more 
pronounced central peak and thinner tails, than would be seen in a normal distribution (Figure 2, Table 1). 

I believe these patterns of skewness and kurtosis are informative in relation to the selective 
pressures on arrival date. The piling up of values on the left side of the distribution implies that the 
pressure to delay migration until food resources are available is stronger than the pressure to arrive 
early to procure a favorable territory. The leptokurtic distribution suggests that the median arrival dates 
are even more advantageous than one would expect if the data were normally distributed. 

Foraging type of birds had an impact on the skewness of yearly arrival date distributions (Table 2). 
I expected that species whose food is more seasonal (herbivorous insects for leaf-gleaners, nectar for 
nectarivores) would show stronger right skew than raptors. That pattern was observed. I also found a 
significant relationship between foraging type and kurtosis. Highest kurtosis values were found for 
nectarivores, suggesting that there is a narrow optimal arrival window for these species. 

The month of arrival was correlated with skewness and kurtosis values. Arrivals in March are more 
perilous than later arrivals because of the fickle nature of the Maine spring. One predicts more of a 
right-skewed and leptokurtic distribution for early arrivals. That pattern was observed with the two 
earliest migrants. Red-winged Blackbird and Common Grackle, showing highest skewness and 
kurtosis [27] presented data to show that long-distance migrants are showing a stronger effect to global 
climate change than short-distance migrants (although see [19]). In New York and Massachusetts, 
Neotropical migrants are arriving about 13 days earlier on average now compared to 50 years ago 
while short-distance migrants are arriving only four days earlier. Insofar as a wintering bird in South 
America cannot predict the weather in the northern United States, we might expect such species to 
show less pronounced skewness and kurtosis than intracontinental migrants that can time their 
migration to regional temperature regimes. Surprisingly, wintering area did not have any statistical 
effect on either skewness or kurtosis. 

For all species, the skewness of yearly distributions of arrival dates was not related to the deviation 
of springtime temperatures from the long-time average. Temperature patterns occur on a broad-scale 
from Delaware north through Maine [19], allowing northward migrants to predict springtime conditions 
500 km or more to the north. One might expect migratory birds to exercise caution when springs are 
colder than normal. I am surprised that none of the 65 species showed this type of response. Similarly, 
none of the species showed a response of yearly kurtosis values to deviation from normal temperatures. 

The larger scale driver of weather, the North Atlantic Oscillation, is related to the timing of bird 
migration [19,28,29]. However, the NAO had no statistically significant relationship with skewness 
and kurtosis for any of the 65 species. 

We have abundant evidence now that both photoperiod and weather patterns influence the timing of 
avian migration [30-33]. Breaking down the distributions of arrival dates for each year jdelds insight 
into the strength of these two cues. As shown in Table 2, stronger correlations consistently exist for 
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adjacent portions of the distribution (e.g., the 5% median with the 10% median) than for portions of 
the distribution more widely separated (e.g., the 10% median with the 75%) median). I interpret this 
consistent result as evidence for the primacy of weather -temperature responses to the migration. If 
photoperiod were driving the migration, we would expect more consistent correlations among portions 
of the distribution more distantly removed. 

Figure 3 explores the relationship of arrival dates between species that arrive during the same time 
period (March-April, April, April-May, May). Because many of the species arriving in the same 
month have similar dietary requirements, correlations of arrival dates are not surprising. The 
unexpected result is the consistently higher proportion of significant correlations between species for 
earlier migrant. The proportion of significant correlations is greatest for the 25% median and least for 
the 75%) median. This result is consistent with the general pattern (Table 1) of a right skew of the 
distributions of arrival dates. The correlation data and the skewness data suggest that a greater penalty 
ensues for arriving too early than for arriving too late. 

5. Conclusions 

In this paper, I have taken a longer, exploratory view of an 18 -year dataset of arrival dates for 
65 species of migratory birds that nest in Maine. Right-skewed arrival dates (Table 1) indicate that 
there is greater selection pressure against arriving too early than arriving later. Leptokurtic patterns 
suggest strong selection for the optimal arrival date. Tighter correlations between the arrival of species 
at the vanguard of migration (the 25% of median) compared to later groups of arrivals supports the 
notion that selection is greatest for the earliest arrivals. Breaking the distribution of arrivals into 
different medians indicates that fine-tuning of the wave of migration occurs as individuals of a species 
migrate into Maine (Table 2). 
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